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NATIONAL FOREWORD 

This Indian Standard (Part 2) which is identical with ISO 6974-2 : 2001 'Natural gas — Determination 
of composition with defined uncertainty by gas chromatography — Part 2: Measuring-system 
characteristics and statistics for processing of data' issued by the International Organization for 
Standardization (ISO) was adopted by the Bureau of Indian Standards on the recommendation of the 
Petroleum Lubricants and Their Related Products Sectional Committee and approval of the 
Petroleum, Coal and Related Products Division Council. 

The text of ISO Standard has been approved as suitable for publication as an Indian Standard without 
deviations. Certain conventions are, however, not identical to those used in Indian Standards. 
Attention is particularly drawn to the following: 

a) Wherever the words 'International Standard' appear referring to this standard, they should 
be read as 'Indian Standard'. 

b) Comma (,) has been used as a decimal marker, while in Indian Standards, the current 
practice is to use a point (.) as the decimal marker. 

In this adopted standard, reference appears to certain International Standards for which Indian 
Standards also exist. The corresponding Indian Standards which are to be substituted in their places 
are listed below along with their degree of equivalence for the editions indicated: 



international Standard 

ISO 6974-1 : 2000 Natural gas — 
Determination of composition with 
defined uncertainty by gas 
chromatography — Part 1: Guidelines 
for tailored analysis 

ISO 6976 : 1995 Natural gas — 
Calculation of calorific values, density, 
relative density and Wobbe index from 
composition 



Corresponding Indian Standard 

IS 15130 (Part 1) : 2002 Natural gas — 
Determination of composition with defined 
uncertainty by gas chromatography: Part 1 
Guidelines for tailored analysis 

IS 14504 : 1998 Natural gas — Calculation 
of calorific values, density, relative density 
and Wobbe index from composition 



Degree of 
Equivalence 

Identical 



do 



For the purpose of deciding whether a particular requirement of this standard is complied with, the 
final value, observed or calculated, expressing the result of a test or analysis, shall be rounded off in 
accordance with IS 2 : 1960 'Rules for rounding off numerical values (revised)'. The number of 
significant places retained in the rounded off value should be same as that of the specified value in 
this standard. 
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Indian Standard 

NATURAL GAS — DETERMINATION OF 

COMPOSITION WITH DEFINED UNCERTAINTY BY 

GAS CHROMATOGRAPHY 

PART 2 MEASURING-SYSTEM CHARACTERISTICS AND STATISTICS FOR 

PROCESSING OF DATA 



1 Scope 

This part of ISO 6974 describes the data processing for the tailored analysis of natural gas. It includes the 
determination of. the measuring system characteristics and the statistical approach to data handling and error 
calculation with the aim of defining the uncertainty in the mole fractions of the component measured. 

This part of ISO 6974 is only applicable in conjunction with part 1 of ISO 6974. 

2 Normative references 

The following normative documents contain provisions which, through reference in this text, constitute provisions of 
this part of ISO 6974. For dated references, subsequent amendments to, or revisions of, any of these publications 
do not apply. However, parties to agreements based on this part of ISO 6974 are encouraged to investigate the 
possibility of applying the most recent editions of the normative documents indicated below. For undated 
references, the latest edition of the normative document referred to applies. Members of ISO and IEC maintain 
registers of currently valid International Standards. 

ISO 6974-1, Natural gas— Determination of composition with. defined uncertainty by gas chromatography — 
Part 1: Guidelines for tailored analysis. 

ISO 6976, Natural gas — Calculation of calorific values, density, relative density and Wobbe index from 

composition. 

3 Terms and definitions 

For the purposes of this part of ISO 6974, the following terms and definitions apply. 

3.1 
response 

output signal of the measuring system for a component that is measured as peak area or peak height 

NOTE Response is expressed in counts. 

3.2 
uncertainty 

estimate attached to a measurement result which characterizes the range of values within which the true value is 
asserted to lie 

NOTE In general, the uncertainty of measurements comprises many components. Some of these components may be 

estimated on the basis of statistical distribution of the results of series of measurements and can be characterized by 
experimental standard deviation. The estimates of other components can only be based on experience or other information. 
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3.3 

certified-reference gas mixture 

CRM 

mixture which is used for the determination of the response curves of the measuring system 

NOTE Certified-reference gas mixtures may be prepared gravimetrically in accordance with ISO 6142l 1 l or ISO 13275I 2 ] or 

certified and validated by comparison with primary standard gas mixtures of closely related composition in accordance with 
ISO 6143P1 (see ISO 141 1 1K]). 

3.4 

working-reference gas mixture 

WRM 

mixture which is used as a working standard for regular calibration of the measuring system 

NOTE Working-reference gas mixtures may be prepared by a gravimetric method in accordance with IS0 6142I 1 ! or 

certified and validated by comparison with CRM of closely related composition in accordance with ISO 6143f 3 ). 

3.5 

direct measurement 

measurement by which individual components and/or groups of components are determined by comparison with 
identical components in the working-reference gas mixture 

3.6 

indirect measurement 

measurement by which individual components and/or groups of components which are themselves not present in 
the working-reference gas mixture are determined using relative response factors to a reference component in the 
working-reference gas mixture 



4 Symbols and subscripts 

4.1 Symbols 

a, /?, < , d, polynomial coefficients of the component /' 

« mC] b mc etc. polynomial coefficients of the components or groups of components to be analysed using direct 
measurement as a function of identical components or groups in the calibration gas 

/i s number of replicates of the sample 

/; wrm number of replicates of the working-reference gas mixture 

K relative response factor with respect to the reference component 

m order of the response function fitted 

n total number of analyses used in the regression 

q number of direct components plus number of indirect components 

R instrument response, expressed as raw data, e.g. area counts or peak height 

flf response factor (based on a single point calibration) 

i MSE residual mean square (MSE) 

5 MS r mean square due to regression (MSR) 
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j§se residual sum of squares (SSE) 

5§sr sum of squares due to regression (SSR) 

s(x) standard deviation of the predicted value 

s(x) standard deviation of a normalized mole fraction 

s{x*) standard deviation of a non-normalized mole fraction 

*(*c,wrc) standard deviation of a mole fraction in the WRM (from the certificate of WRM) 

s(R) standard deviation of the responses at one mole fraction level 

T, the difference in slope between the tangent of the optimal response function in R wi and calibration 

line based on the single point method of the WRM only 

t reduced variable of the r-distribution 

*cwrm mo|G fraction of the WRM (given on the certificate) 

x predicted value of the mole fraction according to a response function 

x* non-normalized mole fraction 

v normalized mole fraction 

.VL,wr lowest mole fraction in the working range 

■»u,wr highest mole fraction in the working range 

i' number of degrees of freedom 

4.2 Subscripts 



c 


certified 


i 


component i 


i 


mole fraction level j 


h 


number of replicates h 


L 


lowest 



mc main components or groups of components to be analysed using direct measurement as a function 

of identical components or groups in the calibration gas 

MSE residual mean square 

MSR mean square due to regression 

oc other components that are not measured and/or can be regarded as being present at a constant 

mole fraction 
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rrf components or groups of components to be analysed using indirect measurement as a function of 

a different reference component in the calibration gas 



ref 


reference component 


SSE 


residual sum of squares 


SSR 


sum of squares due to regression 


s 


sample 


U 


highest (upper limit) 


wr 


working range 


wrm 


working-reference gas mixture 



5 Procedure 

5.1 Step 1: Determination of the optimal response functions 

5.1.1 General remarks 

Detector response functions are determined using certified reference gas mixtures (CRMs). 

For each of the direct components, select the most optimal detector response function which gives the best 
relationship between response signal and component quantity of the CRM. Carry out this step in the following 
cases: 

— initial installation of the system by the supplier; 

— return to operation after a maintenance of the system; 

— return to operation after replacement of a major part of the system, e.g. injection valve, column, detector; 

— return to operation after failure to pass a routine check of the system using test gas(es). 

NOTE , It may be necessary to determine the best statistical curve fit of the results. Although quadratic and cubic curves 
may give a better fit to the data, for a series of gases this may not be feasible. It is more important to obtain consistent results 
for the different gases. In practice, this consistency can often be obtained by simply drawing the relative-response 
(= absolute response/mole fraction) versus mole-fraction curves. 

5.1 .2 Selection of reference gases 

To define the response functions, different mole fraction levels of the direct components are analysed. The number 
of levels needed depends on the order of the response curve that has to be fitted. The minimum number of levels 
needed is the same as the number of polynomial coefficients (see below) in the equation to be fitted. However it is 
recommended to use a few more mole fraction levels than necessary, in most cases, the order of the response 
curve is unknown beforehand but by increasing the number of data points, the response curve is less influenced by 
measurement error. Therefore, a minimum of seven levels is recommended although another number of mole 
fraction levels may be sufficient to obtain a good curve fit. 

Choose these levels so as to spread them uniformly across the specified working range of each component and to 
cover a range which is slightly larger than the anticipated mole fraction range. This can be made possible by using a 
set of multi-component mixtures, each containing different mole fractions of all the direct components. 

NOTE 1 The components can possibly be tested as binary mixtures in methane. However, this would greatly lengthen the 
testing process, as a separate set of cylinders would be needed for each non-methane component. 
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NOTE 2 It is obviously not possible to make a mixture with every component at the highest level required, or at the lowest, so 
most of the multi-component reference gas mixtures will have compositions which are quite different from normal natural gases. 
Provided that the mixtures are shown to be stable in storage and use, this will not create a problem. 

NOTE 3 Prior to the use of a WRM it is necessary to prove that there is no contradiction between the multi-point calibration, 
made by CRMs, and the mole fraction values of the WRM. 

5.1.3 Measurement procedure 

Analyse each of the CRMs, containing the direct components, at least twice (preferably 5 to 10 times) using the gas 
chromatograph. Tabulate the data for the individual responses R;j h of each component i, at each mole fraction; and 
for each replicate h. Before proceeding, inspect the data for outliers. Inspect the results of each group by 
component and level from one calibration interval for outliers or stragglers using Grubb's test or some other 
suitable outlier test. Reject the possible determined outliers. 

NOTE In any set of data, individual results may be found which are not consistent with the other data of that set. These are 

regarded as outliers or stragglers and shall, under the correct circumstances, be eliminated from the set of data. Possible 
causes of outliers can be identified, for example transcription errors or carryover from the previous CRM which was not fully 
purged from the system before the results from the new CRM were recorded. 

5.1 .4 Selectipn of the response function 

5.1.4.1 Introduction 

The response function shall be selected from the following polynomial functions. If a more complex function is 
needed, this means that the method is unsuitable for this application. 

jr, = (a l +)b l R i first order (1) 

Xj = (a, +) bjRj + c t Rf second order (2) 

.*,- = (a,- +) biRj + c t Rf + d,R? third order (3) 

The a-, term (intercept) is placed between brackets because the intercept is tested independently of the order of the 
model, so each of the above listed response functions may also appear without the intercept term. 

Choose the response function with the following assumptions in mind: 

— the most likely response is a straight line through the origin; 

— a zero intercept, implying no response to a component in its absence, is likely for most components; 

— curvature of the response function can be described by a second order term possibly by a third order term; 

— it is unlikely that a higher order term will have greater significance than a lower order term in the same equation; 

— the presence of a maximum or a minimum in the plot of calculated mole fractions versus response within the 
working range is not acceptable. 

5.1.4.2 Regression analysis 

For the data obtained in 5.1 .3, carry out regression analysis using a computer program capable of: 

— performing a least squares regression of component mole fractions, x it as the dependent variable as a function 
of the individual responses, R,j h , as the independent variable; 

— performing a least squares regression of jc,- as a function of R^ and /?|,, ; 
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— performing a least squares regression of x t as a function of R ijh , Rf jh and Rf jh ; 

— producing output for each of the fitted response functions of the: 

— confidence interval for each coefficients; 

— sum of squares due to regression (SSR), j§sr > 

— residual sum of squares (SSE), s|se > 

— mean square due to regression (MSR), %isr! 

— residual mean square (MSE), smse! 

— number of degrees of freedom (v) associated with s^sr ar) d %ise: 

— predicted values of the component mole fractions (Jc, ) ; 

— standard deviations of the predicted values [•*(*,■)] • 

5.1.4.3 Calculation 

Calculate the significance of each of the response functions for each of the components using the following 
procedure (see reference [5] in the Bibliography). 

Calculate the value of t (critical values of t are given in Table A.1 of annex A) for the first order regression /(1), 
relating to b h as 



V 5 mse0) 
The order of the regression, given in parentheses is as follows: 

(1 ) refers to data from the first order regression; 

(2) refers to data from the second order regression; 

(3) refers to data from the third order regression; 

(4) refers to data from the fourth order regression. 

From the first and second order data, calculate r(2), relating to c„ as 



, (2) J 4sp(2)-*Js R< 1) 

V *mse(2) ' 

From the second and third order data, calculate /(3), relating to d h as 



f ( 3 )-J 7^. 6) 

V 5 mse( 3 ) v ' 
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In the acceptance of the instrument carry out a test for the determination of the third order data using a significance 
test between third and fourth order: 



t U\ = HsR(4)-*SSR(3) 
U "V *MSE<4) 



If the fourth order is significant, the system is unsuitable for the application. 

NOTE The uncertainty of the relative, response factor K, compared to the uncertainty of the reference component, is 

considered negligible. [See equations (13) and (15}](A", compared to R rrf ,). 

5.1.4.4 Comparison 

Compare each calculated value of t with the critical values (annex A) at its associated number of degrees of 
freedom. Calculate the number of degrees of freedom using the following equation: 

v= n-(m + 1) 

NOTE For example with seven CRMs each analysed 5 times, n - 5 x 7 = 35, for r(1) v = 35 - (1 + 1) = 33, for /(2) 

K=35-(2+ 1) =32, for/(3) v =31. 

If the calculated value of t is less than the critical value (with the same v), the coefficient is not significant (at a 
significance level of 95 %). If the calculated value of / is greater than the critical value, the coefficient is significant 
(at a significance level of 95 %). 

Select the most appropriate order which includes an intercept term as follows: 

— if r(3) is significant, select the third order response function; 

— if ?(3) is not significant and r(2) is significant, select the second order response function; 

— if /(3) and /(2) are not significant and r(1) is significant, select the first order response function; 

— if /(3), r(2) and /(1) are not significant, do not continue with this procedure as there is no suitable relationship 
between the mole fraction and the response. 

At this stage, consider the significance of the intercept term. From the selected response function, check the 95 % 
confidence interval of the intercept. If this interval includes zero, then the intercept can be set to zero and removed 
from the response function. Repeat the regression analysis for the selected order and for any lower orders. Use the 
same procedure as described above without the intercept (a, - 0). The coefficients calculated by this process are 

identified as b 0l , c 0t and d 0l . 

Calculate r (1), relating to b Qi , with the appropriate .v| SR0 (sum of squares due to zero intercept regression) and 
v mse,o (residual mean square due to zero intercept), using the following equations: 



,(1) = P|1| (7) 

Pmse.oO) 

If a higher order response function has been selected with a zero intercept, then use the zero intercept data to 
calculate t (2) [and /rj( 3 )]- Compare each calculated value of / with the critical values (annex A) at its associated 
number of degrees of freedom. The number of degrees of freedom is calculated as follows: 

v = n — m 

Select the most appropriate order with zero intercept term as follows: 

— if r (3) is significant, select the third order response function; 

— if t {3) is not significant and r (2) is significant, select the second order response function; 
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— if ?o(3) and r (2) are not significant and /o(1) is significant, select the first order response function; 

— if r (3), f (2) and r 0) are not significant, do not continue with this procedure as there is no suitable relationship 
between the mole fraction and the response. 

5.2 Step 2: Selection of the response functions to be used in practice 

5.2.1 General remarks 

The response functions selected in step 1 (see 5.1) are the optimal relationships between mole fraction and 
response. However, it may be necessary to select another response function when the computer program of the 
gas chromatograph is not capable of handling third order functions or functions with an intercept. In this case, a first 
order equation without an intercept has to be used. This approach is usually called the single point calibration 
method. When there is a difference between the optimal function and the first order equation without an intercept 
function, it is obvious that additional uncertainty may be introduced. 

The two methods, namely method A using the optimal function and method B using single point calibration, are 
different in terms of calculation of the uncertainty and in the way calibrations are dealt with by the gas 
chromatograph. Annex B describes how data from a typical test result are handled using method A and method B. 

5.2.2 Method A 

Calculate the desired response functions using CRMs and the coefficients of the response functions. The routine 
calibration, using a WRM, updates the response function coefficients. Possible additional uncertainty (resulting 
from other response functions than the optimal ones) shall be taken into account in the standard calculations. It is 
sufficient to calculate the response functions which are needed in practice. 

5.2.3 Method B 

5.2.3.1 General assumptions 

Calibrate the gas chromatograph using a WRM. This implies the use of a linear calibration (first order regression) 
without an intercept. This procedure does not update a calibration based on CRMs but stands on its own; an 
additional calculation step is needed to take into account the possible additional uncertainty (s wm<i ). 

For this calculation the following assumptions are made: 

— the optimal response functions have already been calculated using step 1 (see 5.1); 

— the mole fractions of the components in the WRM are within the working range of the components in the 
process gas (centre of the range, if possible); 



the response on the WRM is known {R, 



wrm,/7> 



— the working ranges of the components in the process gas are within the application range of the respective 
part of ISO 6974; 

— the mole fractions of the components in the process gas have a more or less normal distribution with the mole 
fraction of the WRM as the average. 

5.2.3.2 Calculations 

Calculate the derivative of the optimal response functions towards the response of the gas chromatograph (/?,-) as 

follows: 

a) optimal response function determined by method A is: 
x, = a, + bjRj + c,Rf + d,Rf 
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b) using single point calibration, differentiate the optimal response function using the following equation: 
x' t = b, + 2c,Rj + 3d,R? 

c) calculate x', at /? wrrTli , (average of two analyses at least) using the following equation: 
x;(wrm) = b t + 2c,-/? wrrTV + 3diRl mi 

d) calculate the response factor based on the WRM using the following equation: 

"•f,wrm,i 



1 wrm.i 



R, 



(8) 



e) calculate 7, by subtracting /?f iWr m,i from x,'(wrm) using the following equation: 
T t = .v,'(wrm) - *f,wrm,/ 



(9) 



NOTE 7, is the difference in the slope between the tangent of the optimal response function in /? W rm,/ ar, d calibration line 

based on the single point method of the WRM only. The working ranges of the components in the process gas should be known 
so as to calculate the extra uncertainty. Divide the working ranges by 4 to calculate the component standard deviations, as 

follows: 



A 'U,wr,, 



L,wr,/ 



(10) 

Moreover multiply s m>l by T, to calculate the standard deviation caused by using the WRM calibration line conforming with 

method B as follows: 



■ v B,i - T, ■ A' W r j( " 
NOTE sqj is used in the calculation of the component mole fraction uncertainties in step 5 (see 5.5). 



(11) 



Caution — When step 1 consists only of replicates of one CRM at mole fraction levels (nearly) similar to those of 
the sample to be analysed, then the optimal response functions (step 1) are not really known and the extra 
calculations for method B cannot be carried out {s Bj = 0, by definition). In this case, no estimate can be made using 
a (possible) calibration line conforming to method B and there is a risk of making systematic errors without 
knowing. This risk shall, in this case, be decreased by limiting the deviation between the composition of the 
calibration gas and the composition of the sample as shown in Table 1. In step 5, sb,,- is then zero. 



Table 1 — Allowed tolerances 


Sample 


Calibration gas mixture 


component mole fraction 


deviation of component mole fraction 


% 


% relative to sample mole fraction 


0,001 to 0,1 


±100 


0,1 to 1 


±50 


1 to 10 


±10 


10 to 50 


±5 


50 to 100 


±3 
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5.3 Step 3: Routine calibration with a working-reference gas mixture 

Carry out routine calibration (daily, weekly, etc.) of the gas chromatograph normally using a working-reference gas 
mixture (WRM). The composition of the WRM shall be within the working range of each component and shall 
contain all the direct components. It is recommended to perform the calibration with the WRM at least in duplicate. 

5.4 Step 4: Calculation of non-normalized component mole fractions 

5.4.1 Method A 

Calculate the non-normalized mole fraction of the direct components in the sample as follows: 

* -*c,wrm,/ . MO\ 
■*mc.i-~ x *mc,/ *'*' 

Xvtrmj 

Calculate the non-normalized mole fraction of the indirect components in the sample as follows: 

* Rrri.i * ., «\rf,f , -*c,wrm,ref „ - /-io\ 
xrrtj = Ki X — x Jrmcref = Ki x ~ — x ~ x *™c,ref ( 13 ) 

ffmc.ref rCmc.ref *wrm,re1 

NOTE Relative response factors, K, for TCD and FID detectors are given for information in annex C. 

5.4.2 Method B 

Calculate the non-normalized mole fraction of the direct components in the sample as follows: 



* _ A c,wrm,i / 14 n 

■Vmc.i ~ A a mc,f \ ' ^> 

Calculate the non-normalized mole fraction of the indirect components in the sample as follows: 

jrrrf ./= K, x X x mc ,ref = Ki x X /? rrf , (15) 

f?mc,ret /vwrm.ref 

NOTE Relative response factors, K, for TCD and FID detectors are given for information in annex C. 

5.5 Step 5: Calculation of non-normalized component standard deviations 

5.5.1 General 

Calculate the standard deviation of the component mole fraction for each of the direct and indirect components. 

5.5.2 Direct components 
5.5.2.1 Standard deviation 

Calculate the standard deviation of the non-normalized mole fraction of the direct components in the sample by 
differentiating equation (12) or (14). This method is the standard method for error propagation. Detailed information 
about this method can be found in literature on this subject (see reference [5] in the Bibliography). 

NOTE In equation (12) there are three sources of error for the non-normalized mole fraction of the direct component (x' mc ). 

Normally, for two of these three sources of error, the contribution of uncertainty to the measurement of the mole fraction of a 
component in the sample is calculated. It is assumed that the contribution of the uncertainty to the mole fraction measurement of 
the components in the WRM [s(jt wrm )] is negligible. If this contribution is not negligible (for example, if the mole fractions in the 
WRM are based on analysis instead of gravimetrical preparation), then equations (17) or (19) should preferably be used to 
calculate this uncertainty. 
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5.5.2.2 Uncertainty determination for method A 

Calculate the uncertainty based on the standard method for error propagation of the non-normalized direct 
components as follows: 



s \Xmc,i>~ x mc,i X ' 



5 \*mc,i) 



Jtmc.i 



2 r 



■H-^wrm,/) 



•*wrm,z 



(16) 



In equation (16) various values determined during the regression (step 1 or step 2) and during the calibration with 
the CRM (step 3) are necessary. Calculate the uncertainty based on the mole fractions of the components in the 
WRM using the following equation: 



s \Xmc,i) ~ -*mc,/ x • 



s \Xmc,i) 



-*mc,; 



s(x 



c,wrm,i i 



l c,wrrru 



(17) 



5.5.2.3 Uncertainty determination for method B 

Calculate the uncertainty based on the standard method for error propagation of the non-normalized direct 
components as follows: 



s{ 



■Xmc.i) — 



s MSE,i x (hwm + hs) 



<wrm 



x/i« 



(18) 



The s MSE value, used in equation (18), is derived from step 1, and is used as an overall approximation of the 
measurement error. The number of replicates of the sample (h s ) is, for example in the case of a laboratory gas 
chromatograph, two or three and in case of a field gas chromatograph, one (by definition). Calculate the uncertainty 
based on the mole fractions of the components in the WRM using the following equation: 



-*( 



xmc,il~ xrnc.i 



s (xmc,i) 



xmc.i 



2 r 



■ *v*x,wrm,; J 

-^cwrm,! 



(19) 



Calculate the uncertainty introduced by using a single point calibration method using the following equation: 



/ * \ it* \2 2 
s Umc,/)=-y/ s Urnc,;) + S Bi 



(20) 



5.5.3 Indirect components 



5.5,3.1 General 



The calculation of the mole fractions of the indirect components involves the use of a reference component (or 
bridge component). As can be seen in equation (13) or (15), this involves the use of extra sources of error in the 

calculation. 



5.5.3.2 Method A 

For method A, account the uncertainty for the sources of error of the indirect components in the sample as follows: 



■ s Urrf,/)-.*rrf,i X ' 



^(^mc.ref) 



-*mc,ref 



J (.*WRM,ref) 
-*WRM,ref 


2 

+ 


*(rtrrf,i) 

KrrU 


2 

+ 



£V_^mc i ref) 
^mc.ref 



(21) 
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Take into account the additional uncertainty based on the mole fractions of the components in the WRM as follows 
s (x T rt,i) = xi-rf ,, x 



s (x*tt,i) 



Xrrli 



2 r 



s \ x c,vjrw,i ) 



-^wrm.i 



(22) 



5.5.3.3 Method B 



For method B, account the uncertainty for the sources of error in the mole fractions of the indirect components in 
the sample as follows: 



, * v j- s MSE,ref x C*wrm + M 

s Urrf.i) = J 7 —7 

V "wrm x "s 



(23) 



Derive the s MSE value, used in equation (23), from step 1, and use this as an overall approximation of the 
measurement error. The number of replicates of the sample (h s ) is, for example in the case of a laboratory gas 
chromatograph, two or three and in case of a field gas chromatograph, one (Dy definition). 



Take into account the additional uncertainty based on the mole fractions of the components in the WRM as follows: 

(24) 



■rUrrf,/) = *rrf,i x - 



s(x*rt,) 



xrri.i 



s(x 



c,wrm,( 



l c,wrm,i 



Take into account the uncertainty introduced by using a single point calibration method as follows: 



s ( x rf,i)=\j s(x*rri,i) + 4,ref 



5.6 Step 6: Calculation of normalized component mole fractions 



(25) 



In the calculation of the normalized mole fractions, the mole fractions of the "other components" (with a fixed mole 
fraction) are also needed. It is assumed that the total amount of these components in the sample is small compared 
to the rest of the (direct and indirect) components. 

Normalization of the results is allowed if the sum of the mole fractions of the non-normalized component mole 
fractions is no less than 0,98 or greater than 1 ,02. Calculate the fraction of the normalized components using the 
following equation: 



*i = 



X( 



I * 

Zj Xw 

w = 1 



x (1 - Jtoc) 



(26) 



5.7 Step 7: Calculation of normalized component standard deviations 

Calculate the normalized component standard deviations of the normalized components using the following 
equations. 



*(xi) = x, x V-%2- x ,(^) 2 + j^s( x lr 

V Xi w=1 



(27) 
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At the right side of equation (27), the sum of alt the squared standard deviations (variances) of the non-normalized 
component mole fractions (direct + indirect components) is calculated. Equation (27) is an approximation and is 
based on the assumption that the sum of the non-normalized mole fractions is close to 1 . 

NOTE The (summed) mole fraction of the other components (x oc ) is, in most cases, very small when compared to the rest 

of the components. Therefore the standard deviation of the other components (jc oc ) is not taken into account in the calculation of 
the standard deviation of the normalized components. This would make the calculations much more complicated without having 
any significant impact on the end result. 

5.8 Step 8: Calculation of uncertainty and repeatability of the results 

5.8.1 Absolute uncertainty 

Approximate the absolute uncertainty, U abs (x,), of the normalized components using the following equations. 

U abs (xi) = txs( Xl ) (28) 

t can be obtained from the Table A.1 in annex A. The number of degrees of freedom of this t is derived from step 1 
[h-(w-f1) in case of a response function with intercept, and n-m in case of a response function without 
intercept]. 

5.8.2 Relative uncertainty 

Calculate the relative uncertainty £/ re i(*,), expressed as a percentage, of the normalized components using the 
following equation: 

f/relU,)- ^ absU) x100 (29) 

Xj 

5.9 Relationship between ISO 6974-2 and ISO 6976 

5.9.1 General 

In ISO 6976, data are needed to calculate the listed physical properties and their repeatability. In this subclause, a 
link is established between the output of this part of ISO 6974 and ISO 6976. 

5.9.2 Non-normalized mole fraction 

Calculate the repeatability of the standard deviations of the non-normalized components mole fractions as follows: 

r* = 2y[2xs{ x *) 
where 

Ajr/ = n 

5.9.3 Normalized mole fraction 

Calculate the repeatability of the normalized components mole fractions as follows: 

r, = 2^xs(x,) ^ 

where 

xj is the normalized mole fraction 

zlv, = r, 
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Annex A 

(informative) 

Critical values of t 



The critical values of t are given for confidence levels of 95 % and 99 % in Table A.1. 



Table A.1 — Critical values of / 



Degrees of freedom 


Confidence level 


95% 


99% 


1 


12,7 


63,7 


2 


4,30 


9,92 


3 


3,18 


5,84 


4 


2,78 


4,60 


5 


2,57 


4,03 


6 


2,45 


3,71 


7 


2,36 


3,50 


8 


2,31 


3,36 


9 


2,26 


3,25 


10 


2,23 


3,17 


11 


2,20 


3,11 


12 


2,18 


3,05 


13 


2,16 


3,01 


14 


2,14 


2,98 


15 


2,13 


2,95 


16 


2,12 


2,92 


17 


2,11 


2,90 


18 


2,10 


2,88 


19 


2,09 


2,86 


20 


2,09 


2,85 


CO 


1,96 


2,58 
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Annex B 

(informative) 

Example 



This annex gives, as an example, an overview of experimental results obtained with a gas chromatograph. Data 
from analyses are presented in Table B.1 and treatment of the results following the steps outlined in clause 5 are 
presented in Tables B.2 to B.10. 

Table B.1 — Raw data from analyses using certified reference gas mixtures (CRM) (triplicate analysis) 





Gas 1 


Gas 2 


Methane, CRM 
Gas 3 Gas 4 


Gas 5 


Gas 6 


Gas 7 


Mole fraction (%) 


65,146 


85,776 


75,602 


88,766 


79,758 


70,301 


95,203 


Responses 


165798,87 
165 979,54 
165981,30 


214424,77 
214555,59 
214285,16 


190392,92 
190208,11 
190246,59 


221 549,75 
221 426,74 
221 133,67 


200286,27 
200304,79 
200070,49 


177903,69 
177911,01 
177783,43 


236194,26 
236139,76 
236314,58 




Gas 1 


Gas 2 


Ethane, C 
Gas 3 


RM 
Gas 4 


Gas 5 


Gas 6 


Gas 7 


Mole fraction (%) 


5,702 


3,439 


10,894 


0,595 


8,374 


7,380 


1,451 


Responses 


23 570,03 
23 585,87 
23 590,53 


14 283,97 
14295,91 
14275,36 


44673,11 
44659,86 
44625,63 


2 502,28 
2 496,92 
2 494,49 


34442,71 
34428,37 
34403,68 


30387,47 
30389,68 
30378,11 


6 044,86 
6044,94 
6048,77 




Gas 1 


Gas 2 


Propane, ( 
Gas 3 


:rm 

Gas 4 


Gas 5 


Gas 6 


Gas 7 


Mole fraction (%) 


3,927 


3,422 


0,081 


2,291 


0,820 


2,872 


0,899 


Responses 


20659,79 
20668,98 
20 680,61 


18044,84 
18064,06 
18041,95 


439,66 
434,65 
434,00 


12152,91 
12154,09 
12132,08 


4269,36 
4267,50 
4266,81 


15084,82 
15133,59 
15081,18 


4791,53 
4795,68 
4 828,73 




Gas1 


2- 
Gas2 


Methylpropa 
Gas 3 


ne, CRM 
Gas 4 


Gas 5 


Gas 6 


Gas 7 


Mole fraction (%) 


0,033 


0,144 


0,230 


0,366 


0,472 


0,592 


0,558 


Responses 


212,41 
214,33 
214,28 


923,47 
925,21 
924,62 


1 460,75 
1454,15 
1 455,57 


2304,16 
2 321,25 
2294,51 


2 958,62 
2944,19 
2 940,71 


3681,85 
3665,56 
3680,96 


3515,96 
3519,30 
3521,76 
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Table B.1 (continued) 





Gas 1 


Gas 2 


n-Butane, CRM 
Gas 3 Gas 4 


Gas 5 


Gas 6 


Gas 7 


Mole fraction (%) 


0,030 


0,143 


0,235 


0,348 


0,460 


0,685 


0,621 


Responses 


198,80 
202,53 
202,37 


873,20 
874,69 
872,91 


1 460,75 
1 453,74 
1456,10 


2187,50 
2160,36 
2181,44 


2842,89 
2876,33 
2872,86 


4298,82 
4273,51 
4298,56 


3826,17 
3828,50 
3828,75 




Gas1 


Gas 2 


Nitrogen, CRM 
Gas 3 Gas 4 


Gas 5 


Gas 6 


Gas 7 


Mole fraction (%) 


17,605 


2,481 


12,733 


5,751 


4,326 


8,852 


0,301 


Responses 


53439,93 
53471,20 
53475,91 


7757,14 
7760,39 
7756,20 


38806,68 
38767,16 
38780,51 


17839,91 
17831,05 
17819,08 


13432,12 
13425,94 
13417,57 


27199,58 
27231,83 
27219,39 


993,44 
992,75 
994,79 




Gas 1 


C 
Gas 2 


arbon dioxid 
Gas 3 


e,CRM 
Gas 4 


Gas 5 


Gas 6 


Gas 7 


Mole fraction (%) 


7,558 


4,595 


0,225 


1,883 


5,791 


9,317 


0,967 


Responses 


27318,70 
27337,69 
27 348,80 


16645,62 
16658,36 
16634,59 


836,95 
834,69 
835,18 


6837,86 
6834,00 
6829,18 


20938,43 
20939,91 
20919,43 


33 587,92 
33598,91 
33586,73 


3 515,53 
3513,75 
3516,45 


Component 


Wor 
N 2 


king-referenc 
C0 2 


;e gas mixtu 
CH 4 


re (duplicate 
C 2 Hg 


analyses) 
C3H8 


iso- C4H1Q 


/T-C4H-1Q 


Mole fraction (%) 


13,703 


1,049 


82,568 


2,099 


0,431 


0,068 


0,082 


Responses 


41 139,33 
41139,42 


3814,33 
3814,36 


205 395,02 
205395,22 


12101,09 
12101,14 


2276,10 
2276,13 


440,22 
440,24 


513,29 
513,31 


Component 


N 2 


Sam 

C0 2 


}le (duplicate 
CH 4 


; analysis) 
C 2 H 6 


CaHa 


/SO-C4HK) 


JT-C4H10 


Responses 


40831,46 
40823,92 


3808,56 
3807,52 


205 856,65 
205 934,98 


11975,91 
11977,43 


2 285,85 
2286,06 


426,39 
426,93 


529,00 
529,01 


Component 


neo-CsH^ 


/SO-C5H12 


n-C^H-\2 


c 6+ 








Responses 


54,74 
54,43 


144,81 
144,87 


140,46 
140,31 


553,32 
557,18 



Step 1 : Regression analysis (see 5.1.4) 

The response function for carbon dioxide has been selected as the example for calculations. Calculations for the 
other components can be made in a similar way. 

A computer program has been used to fit the first, second and third order response functions (including the 
intercept). The results of these fits are given in Table B.2 (mole fractions have been divided by 100 to be expressed 
as a percentage). 
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Table B.2 — Fits to response function 



Order 


2 
^SSR 


■*MSE 


First 

Second 

Third 


0,021492884 
0,021492970 
0,021492985 


7,22887x10-9 
2,84930x10-9 
2,18136x10-9 



Step 2: Selection of the response function (see 5.1.4) 

The t value calculated for each order of polynomial function is given in Table B.3. 

Table B.3 — Criteria for choice of the response function 



Order 


calculated / value 


V 


critical t value 


First 

Second 

Third 


1 724,297 
5,494 
2,622 


19 
18 
17 


2,09 
2,10 
2,11 



Conclusion: The third order term is significant (t(3) > critical / value), the most appropriate model at this stage is 
the third order response function. 

Verification of the correct function chosen 

Test of the intercept of the third order response function: 

Intercept = -7,541 x 10-5 ± 6,343x10-5 

95 % confidence interval = -1,388x10-4 <; -7,541x10-5 < -1,1 98x1 0" 5 

The 95 % confidence interval does not include zero. Therefore, the intercept is significant. 

Conclusion: For carbon dioxide, the third order model with intercept is selected. 

The optimal response function for each of the other components is calculated for the other components in the 
CRMs. 



Table B.4 — The optimal response functions for the other components in the CRMs 


Component 


Order 


Intercept? 


0| 


bi 


Ci 


di 


Methane 


3 


yes 


-4,126x10-1 


9,745x10-6 


-2,783x10-11 


4,670x10-17 


Ethane 


3 


no 


— 


2,382x10-6 


1,968x10-12 


-1,512x10-17 


Propane 


1 


no 


— 


1,897x10-6 


— 


— 


/'so-Butane 


1 


yes 


-3,337x10-5 


1,607x10-6 


— 


— 


n-Butane 


1 


no 


— 


1,607x10-6 


— 


— 


Nitrogen 


3 


no 


— 


3,155x10-6 


4,919x10-12 


-4,377x10-17 


Carbon dioxide 


3 


yes 


-7,541x10-5 


2,775x10-6 


-1,063x10-12 


3,201x10-17 ■ 
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The relative response factor (K) used for the back-flush measured Cg+ component is estimated to be 0,59 in this 
example. 

Method A: The response functions used will be the same as the optimal response functions calculated in step 1. In 
practice, when method A is followed and other response functions are used, the desired response functions have 
to be fitted to the CRM data. The rest of procedure of method A is exactly the same, independent of the type of 
response functions used in practice. 

No extra calculations are needed, the relevant information is obtained from step 1. 

Method B: The response function used will be first order without the intercept, and the extra uncertainty due to 
using a "non-optimal" response function will be calculated. 

Example: Calculation of carbon dioxide 

— Optimal response function determined using method A: 

x„ n = -7,541x10 -5 + 2,775 x 10 ~ 6 /?„ n - 1,063 x10~ 12 /? 2 + 3,201 x10" 17 /?^„ 

— Using single point calibration (method B), differentiate the optimal response function as follows: 
x' CQ = 2,775x10 6 - 2,126 x10^ 12 ^ co + 9,603 x10~ 17 /? 2 o 

— Response of the WRM for carbon dioxide: /?co 2 = 3814,345 (average of the two analyses). Thus: 
jc' co = 2,775x10^ 6 - 2,126 x10" 12 x3 814,345 + 9,603x10~ 17 x(3 814,345) 2 = 2,7683x 10~ 6 

— Response factor based on the WRM: flf, wr ,co 2 = 0,01049/381 4,345 = 2,7501 x10~ 6 
r C0 2 =y wr,C0 2 ( wr > - R f,w,C0 2 = 2,7683x10~ 6 - 2,750 1x10" 6 = 1,8200x10~ 8 

— The working range for carbon dioxide is for example (at a certain location where this particular gas 
chromatograph is located): 

*L,wr,co2 - 0,005 mole fraction 
*u,wr,co 2 = °.° 20 mo[e fraction 
Thus s W r,C0 2 = (*U,wr,C02 " *L,wr.C0 2 ) I 4 = (0,020 - 0,005) / 4 = 0,00375 mole fraction 

— The extra standard deviation (using method B) is therefore: 

swrm,co 2 = 7co 2 x s m,co 2 = 1,8200 x 10-8 * 0,00375 = 6,825 x 10-n mole fraction 

— For the other components in the CRM see Table B.5 



18 



IS 15130 (Part 2): 2007 
ISO 6974-2 : 2001 



Table B.5 — Calculation for other components in the CRM 



Component 


Ti 


mole fraction 


-*U,wr,i 

mole fraction 


s B,i 

mole fraction 


Methane 


1,9790x10-6 


0,80 


0,84 


1,9790x10-8 


Ethane 


6,6905x10-7 


0,01 


0,03 


3,3453x10-9 


Propane 


3,4200x10-9 


0,002 


0,006 


3,4200x10-12 


/'so-Butane 


6,2350 x 10-8 


0,0005 


0,0006 


6,2350x10-12 


n-Butane 


9,4900 x 10 9 


0,0006 


0,001 


9,4900x10-13 


Nitrogen 


-4,761 4 x 10-8 


0,12 


0,14 


-2,3807x10-10 


Carbon dioxide 


1,8200x10-8 


0,005 


0,02 


6,8250x10-11 


NOTE The standard < 


deviations given in this table are to be used in step 5. 





Step 3: (see 5.3) 

Two analyses have been performed using the WRM. 

Step 4: (see 5.4) 

For the calculation of the non-normalized mole fractions, the average value of the responses of the components in 
the sample and the WRM have been taken for the rest of the calculations. 

Method A 

Direct component (C0 2 ): 
tmc,co 2 = 1 ,0478 x 10 2 mole fraction 

v mc,co 2 = 0,01049/0,010495 x 1,0478 x 10-2 = 0,010473 mole fraction or 1,0473 % (mole fraction expressed as 
a percentage) 

Indirect component (neo-CsH-^. C3 as reference): 

-£mc,c 3 H 8 = 1 - 87 ' 1 9x10- 5 % (mole fraction expressed as a percentage) 

xrrf.neo-c 5 H 12 = 0,007753 % (mole fraction expressed as a percentage) 

Method B 

Direct component (C0 2 ): 

jtmc,co 2 = 1 ,04727 x 10-2 mole fraction or 1 ,0473 % (mole fraction expressed as a percentage) 

Indirect component (neo-CsH^i C3 as reference): 

•*rrf,neo-C 5 H 12 = 7,752 1 x 10~ 5 mole fraction or 0,007752 % (mole fraction expressed as a percentage) 
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Table B.6 — Comparison of the non-normalized mole fractions in the sample 
determined using method A and method B 



Component 


* 
Xi 

mole fraction 


Method A 


Method B 


Methane 


0,82781 


0,82769 


Ethane 


0,020772 


0,020774 


Propane 


0,004329 


0,004329 


/so-Butane 


0,0006580 


0,0006590 


n-Butane 


0,0008451 


0,0008451 


Nitrogen 


0,135 97 


0,13599 


Carbon dioxide 


0,010473 


0,010472 


neo-Pentane 


0,00007752 


0,000077 52 


/so-Pentane 


0,00020021 


0,00020021 


n-Pentane 


0,00019406 


0,000194 06 


c 6+ 


0,000 620 33 


0,000 620 33 



Step 5 (see 5.5.) 

For the calculation of the mole fractions, the average value of the responses of the components in the sample anc 
the WRM have been taken, before making the remainder of the calcuiations. 



Equations (18) and (23) have been used to calculate the standard deviations for method B where lu 



'wrm 



Further it is assumed that the additional uncertainty caused by the WRM is negligible. 

Values for the calculations to determine the standard deviations of non-normalized mole fractions obtained using 
methods A and B are given in Table B.7. 
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Table B.7 — Comparison of the calculation of the standard deviations of the 
non-normalized mole fractions in the sample using method A and method B 





six*) 


Component 


mole fraction 


Method A 


Method B 


Methane 


0,0005753 


0,0005157 


Ethane 


0,00003484 


0,00004199 


Propane 


0,00009337 


0,00009320 


/so-Butane 


0,00003322 


0,00002956 


n-Butane 


0,00003584 


0,00003544 


Nitrogen 


0,0001347 


0,0001100 


Carbon dioxide 


0,00005176 


0,00004671 


neo-Pentane 


0,000001 701 


0,00009320 


/so-Pentane 


0,000004319 


0,00009320 


n-Pentane 


0,000004188 


0,00009320 


c 6+ 


0,00001372 


0,00009320 



Step 6 (see 5.6) 

In this example there are no other components, so jcqc = 0. Mole fractions have been normalized and values are 
given in Table B.8. 

Table B.8 — Comparison of the normalized mole fractions in the sample 
determined using method A and method B 



Component 


Xi 

mole fraction 


Method A 


Method B 


Methane 


0,82619 


0,82616 


Ethane 


0,020732 


0,020735 


Propane 


0,0043202 


0,0043206 


/so-Butane 


0,00065671 


0,00065782 


n- Butane 


0,00084344 


0,00084352 


Nitrogen 


0,13571 


0,13574 


Carbon dioxide 


0,010452 


0,010453 


neo-Pentane 


0,000077369 


0,000077377 


/'so-Pentane 


0,00019982 


0,00019984 


n-Pentane 


0,00019368 


0,00019370 


c 6+ 


0,00061912 


0,00061918 



Step 7 (see 5.7) 

Mole fractions have been normalized and values are given in Table B.8. 
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Table B.9 — Comparison of the calculation of the standard deviations of the 
normalized mole fractions in the sample using method A and method B 



Component 


s( Xi ) 
mole fraction 


Method A 


Method B 


Methane 


0,0001804 


0,0002234 


Ethane 


0,00003627 


0,00004271 


Propane 


0,000092 83 


0,00009266 


/so-Butane 


0,00003313 


0,00002949 


n-Butane 


0,00003574 


0,00003534 


Nitrogen 


0,0001410 


0,0001217 


Carbon dioxide 


0,00005150 


0,00004651 


neo-Pentane 


0,000001698 


0,00009302 


/soPentane 


0,000004311 


0,00009301 


n-Pentane 


0,000004181 


0,00009301 


C6+ 


0,00001369 


0,00009297 



Step 8 (see 5.8) 

The number of degrees of freedom for the different components can be calculated for both methods from the 
results of step 1 . Based on the number of degrees of freedom, the critical t value can be obtained from Table A.1 of 
annex A. This t value is then multiplied by the normalized component standard deviations (in this case). 

Table B.10 — Results of calculation of uncertainty and repeatability 









^abs(^) 


kYelt^ 


</) 


Component 






mole fraction 


mole fraction (%) 


V 


/(Table A.1) 


Method A 


Method B 


Method A 


Method B 


Methane 


17 


2,11 


0,00003807 


0,0004714 


0,04608 


0,05706 


Ethane 


18 


2,10 


0,000076017 


0,000089 69 


0,367 4 


.0,4325 


Propane 


20 


2,09 


0,0001940 


0,0001937 


4,491 


4,482 


/so-Butane 


19 


2,09 


0,00006925 


0,00006163 - 


10,54 


9,368 


n-Butane 


20 


2,09 


0,00007470 


0,00007387 


8,856 


8,757 


Nitrogen 


18 


2,10 


0,0002960 


0,0002656 


0,2181 


0,1883 


Carbon dioxide 


17 


2,11 


0,0001087 


0,00009814 


1,034 


0,938 9 


.neo-Pentane 


20 


2,09 


0,000003549 


0,0001944 


4,587 


251,3 


/so-Pentane 


20 


2,09 


0,000009011 


0,0001944 


4,510 


97,27 


n-Pentane 


20 


2,09 


0,000008738 


0,0001944 


4,512 


100,4 


c 6+ 


20 


2,09 


0,00002862 


0,0001943 


4,622 9 


31,38 
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Annex C 

(informative) 

Relative response factors (K) 



C.1 Relative response factors for FID 

For an FID detector, the relative response factors are calculated as the ratio between the carbon number of the 
reference component to the carbon number of the sample component. Values for FID relative response factors for 
reference components, propane and n-butane, are given in Table C.1. 



Tabl 


e C.1 — Relative response factors 


Component 


Response factor relative to 

propane n-butane 


Propane 


1,000 


1,333 


/so-Butane 


0,750 


1,000 


n-Butane 


0,750 


1,000 


Pentanes 


0,600 


0,800 


Hexanes 


0,500 


0,667 


Heptanes 


0,429 


0,571 


Octanes 


0,375 


0,500 


Benzene 


0,500 


0,667 


Cyclohexane 


0,500 


0,667 


Methyl-cyclohexane 


0,429 


0,571 


Toluene 


0,429 


0,571 



These figures are derived from reference [6] of the Bibliography. 



C.2 Typical relative response factors for TCD 

For a TCD detector, the relative response factors are based on practical measurements and can be used as a 
guideline. An alternative method to obtain TCD relative response factors is, for example, to determine them by an 
(extended) analysis using an FID detector and the figures listed above. Some experimental data are given in 
Table C.2. 

Table C.2 — Experimental K data 



Component 


K 
(relative to propane) 


neoPentane 

/so-Pentane 

n-Pentane 


0,75 
0,73 
0,73 
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